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Abstract 

The spinless Salpeter equation may be considered either as a standard approximation 
to the Bethe-Salpeter formahsm, designed for the description of bound states within a 
relativistic quantum field theory, or as the most simple, to a certain extent relativistic 
generalization of the costumary nonrelativistic Schrodinger formalism. Because of the 
presence of the rather difficult-to-handle square-root operator of the relativistic kinetic 
energy in the corresponding Hamiltonian, very frequently the corresponding (discrete) 
spectrum of energy eigenvalues cannot be determined analytically. Therefore, we show 
how to calculate, by some clever choice of basis vectors in the Hilbert space of solutions, 
for the rather large class of power-law potentials, at least (sometimes excellent!) upper 
bounds on these energy eigenvalues, for the lowest-lying levels this even analytically. 
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1 Introduction: The Spinless Salpeter Equation 

One's attitude to the well-known "spinless Salpeter equation" may be 
reflected by either of the following two approaches (or points of view): 

• On the one hand, this spinless Salpeter equation may be regarded 
to represent some standard approximation to the Bethe-Salpeter 
formalism for the description of bound states within a relativistic 
quantum field theory. It may be derived from the Bethe-Salpeter 
equation by two steps: 

1. Eliminate — in full accordance with the spirit of instantaneous 
interactions — any dependence on timelike variables to obtain 
in this way the so-called Salpeter equation []2j]. 

2. Neglect any reference to all the spin degrees of freedom of the 
involved bound-state constituents and restrict your formalism 
exclusively to positive-energy solutions. 

• On the other hand, this spinless Salpeter equation may be viewed 
as one of the most straightforward generalizations of the standard 
nonrelativistic quantum theory towards the reconciliation with all 
the requirements imposed by special relativity. To be precise, this 
generalization consists of incorporating the square-root operator 
of the relativistic expression for the kinetic energy of the involved 
particles. For the particular case of two particles of equal mass m 
and relative momentum p, the kinetic-energy operator T is given 



All the forces operating between the bound-state constituents are 
tacitly assumed to be described by an arbitrary static interaction 
potential V. For the special case of two particles, this interaction 
potential should depend only on the relative coordinate x of these 
particles: V = V{x.). 

In any case, the self-adjoint Hamiltonian H governing the dynamics of 
any quantum system to be described by the spinless Salpeter equation 
will be of the form 



The two-particle spinless Salpeter equation to be investigated here is 
then nothing else but the eigenvalue problem for this Hamiltonian 



by 







H\xk) = Ek\xk) , ^ = 0,1,2,... 
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for Hilbert-space eigenvectors \xk) corresponding to energy eigenvalues 

^ {Xk\H\xk) 
iXklXk) 

For the sake of simplicity, we shall focus our attention to the physically 
most relevant case of central potentials, i. e., potentials which depend 
only on the modulus |x| of the configuration-space relative coordinate: 

V = V{\^\). (3) 

In the above form, the spinless Salpeter equation appears to be a very 
promising candidate for the (semi)relativistic description of hadrons as 
bound states of (constituent) quarks within the framework of potential 
models g, |^, or, at least, the first step in the correct direction |7|]. 

However, the presence of the relativistic kinetic-energy operator ([T|) 
in (|2]) or, to do justice to the spinless Salpeter equation, the nonlocality 
of this operator that is, more precisely, of either the kinetic-energy 
operator T in configuration space or the interaction-potential operator 
V in momentum space, renders difficult to arrive at rigorous analytical 
statements about the corresponding energy spectrum. In view of this, 
numerous attempts to circumvent these problems have been proposed. 
Some very brief account of the history of these attempts may be found, 
for instance, in Ref. These approaches include, among others, the 
development of elaborate numerical approximation methods [Ty, |TT| 
as well as the construction of effective Hamiltonians which, in spite of 
their apparently nonrelativistic form, incorporate relativistic effects by 
sophisticated momentum dependence of the involved parameters [T^ . 
A lot of information on the solutions of the spinless Salpeter equation 
may even be gained by application of a relativistic virial theorem [jT^ 
most easily derived from a rather general "master virial theorem" [jn 

The (from the physical point of view perhaps most interesting) case 
of a Coulomb-type static interaction potential, the so-called relativistic 
Coulomb problem, has been investigated particularly carefully. For the 
corresponding lowest-lying energy eigenvalues, both lower [jTB], [T^ and 
upper [|TB|, |T7|, |TB|, |T^ bounds have been derived and series expansions 
23] in powers of the involved fine structure constant have been given. 



Here, we intend to pave the way for the calculation of upper bounds 
on the energy eigenvalues of the spinless Salpeter equation with rather 
arbitrary interaction potentials. To this end, we apply the famous min- 
max principle — which controls any such attempt — in a particular basis 
of our trial space, characterized by generalized Laguerre polynomials. 
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2 Minimum— Maximum Principle and Rayleigh— 
Ritz Variational Technique 

The derivation of upper bounds on the eigenvalues of some operator H 
makes, of course, only sense for those operators H which are bounded 
from below. Accordingly, let us assume from now on that the arbitrary 
interaction potential (|^) in our semirelativistic Hamiltonian (g) is such 
that this necessary prerequisite holds. For example, for the crucial case 
of a Coulomb-type static interaction potential, the so-called relativistic 
Coulomb problem, the demanded semi-boundedness of the spectrum of 
the Hamiltonian H has been (rigorously) demonstrated by Herbst []T^ . 



The theoretical basis as well as the primary tool for the derivation of 
rigorous upper bounds on the eigenvalues of some self- adjoint operator 
is, beyond doubt, the so-called min-max principle []2T[]. An immediate 
consequence of this min-max principle is the Rayleigh-Ritz technique: 
Let i7 be a semi-bounded self-adjoint operator. Let Ek, /c = 0, 1, 2, . . ., 
denote the eigenvalues of H, ordered according to Eq < Ei < E2 < . . .. 
Let Dd be some d-dimensional subspace of the domain of H and let Ek^ 
/c = 0, 1, . . . , d — 1, denote all d eigenvalues of this operator H restricted 
to the space D^, ordered according to Eq < Ei < . . . < Ed-i- Then the 
kth eigenvalue Ek (counting multiplicity!^) of H satisfies the inequality 

Ek<Ek, /c = 0,l,...,d-l . 

(For a discussion of the history of inequalities and variational methods 
for eigenvalue problems, see, e. g., Ref. for some applications, see, 
e. g., Ref. 0.) 



Now, let us assume that this d-dimensional subspace Dd is spanned 
by some set of d orthonormalized (and therefore beyond doubt linearly 
independent) basis vectors \i/Jk): = 0, 1, . . . , (i — 1: 

{Hi^j)=^ij: = 0,l,...,(i- 1 . 

Then the set of eigenvalues E may immediately be determined as the d 
roots of the characteristic equation 

det {{HH\ijj) - E6,j) = , J = 0, 1, . . . , d - 1 , (4) 

as becomes clear from an expansion of any eigenvector of the restricted 
operator H in terms of the set of basis vectors \i/Jk), /c = 0, 1, . . . , (i — 1, 
of the subspace Dd- 

^ For instance, for a Hamiltonian H depending only on the moduli of momentum p 
and coordinate x, respectively, states of given orbital angular momentum but different 
projections of the latter will be degenerate. 
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3 Generalized Laguerre Basis 

The crucial step in any investigation of the present type is the suitable 
choice of a basis in the subspace Dd- For the case of the semirelativistic 
Hamiltonian (0), we find it convenient to work in a basis which involves 
the so-called generalized Laguerre polynomials. The latter are specific 
orthogonal polynomials, defined by the power series []2^ 



= t i-K 

r=0 

and normalized according to 




^• 

Consequently, introducing two variational parameters, namely, one, /i, 
with the dimension of mass as well as a dimensionless one, /3, a generic 
trial vector \^JJ) of the subspace D^, with orbital angular momentum i 
and its projection m, will be characterized by the following admittedly 
very suggestive ansatz for its coordinate-space representation i/jk/m{'^)- 

^M-(x) =AA|x|^+^-iexp(-^|x|)4^^(2/i|x|)3^£^(l^x) , (5) 

where normalizability restricts the variational parameter ji to positive 
values, 

/i > . 

Here, yemi^) are the spherical harmonics for angular momentum i and 
projection m depending on the solid angle Q; they are orthonormalized 
according to 

J dQ ^I^(O) ytm'i^) = See Smm' ■ (6) 

The proper orthonormalization of the ansatz (|5]) fixes the parameter 7 
necessarily to the value 7 = 2 £ + 2 and determines the normalization 
constant J\f: 



^r(2£ + 2/3 + A; + l 

X Lf+^^)(2/i|x|)3^,^(n.) 
satisfies the normalization condition 



x| "^'^ exp(— /i|x|) 
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Rather obviously, normahzabihty constrains the variational parameter 
P too, namely, to a range characterized by 2 > — 1, i. e., to the range 

^ 2 

The Fourier transform ipk,imip) of the above trial function involves the 
hypergeometric series F, defined with the help of the gamma function 

r by 



^^^^^^^^ r(w; + n) n! 



it reads 



^r(2£ + 2/3 + /c + i) 2^+1/2 r(^ + 1) 



ro r! I k-r I (p2 ^ ^2)(2£+/3+r+2)/2 



X ^-^;^+2^^2T^J^-(^p) 

and satisfies the normalization condition 

In principle, it is straightforward to calculate the expectation values 



of the Hamiltonian (0), necessary for applying the min-max principle. 
Due to the orthonormalization (g) of the spherical harmonics y£jn{Q), 
however, only matrix elements taken between states of identical orbital 
angular momentum i and its projection m will be nonvanishing. 



4 Power-Law Potentials 

When speculating about the possible shape of a physically meaningful 
(or phenomenologically acceptable) interaction potential, the very first 
idea which unavoidably comes to one's mind as a reasonable candidate 
is an interaction potential of the power-law form, the power being only 
constrained by requiring that the Hamiltonian is bounded from below: 

y(|x|) = E«Jx^ , (7) 
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with sets of arbitrary real constants a„ and the latter only subject 
to the constraint 

&n ^ ~1 if On < . 
By close inspection of our ansatz (|5|) it should become clear that we are 
able to handle even potentials of the type "power-times-exponential," 
that is, potentials of the form 

^(|x|) = E^n |x|^"exp(c^ |x|) , bn>-l if a„ < . 

n 

It is a rather simple task to write down the matrix elements for the 
power-law potential (|7|): 



= J2an [ d^x?/'*^^(x) |x|^"^_^-£^(x) 



^ r(2 ^ + 2 + z + 1) r(2 £ + 2 /3 + J + 1 



E E 



■1)' 



[2iiy-r'^o^o Hs! I i-r 



j + 2^ + 2/3 
J - s 



X r{2£ + 2/3 + bn + r + s + l) . 

For instance, considering merely radial excitations by letting i = and 
choosing, just for the sake of definiteness, for the variational parameter 
(3 the value /3 = I, the explicit form of the potential matrix V = (Vij) 
is 

/ 3 -VSbn 



D n 



r(3 + b„) 



n 

V36„ 3 + 6„ + 62 



5 Analytically Evaluable Special Cases 

It should be really no great surprise that the evaluation of the matrix 
elements of the kinetic-energy operator T, 

is somewhat more delicate than the previous calculation of the matrix 
elements of the power-law potentials V. Consequently, let us focus our 
attention to those situations which allow for a fully analytic evaluation 
of the above kinetic-energy matrix elements. 



7 



5.1 Orbital Excitations 

On the one hand, we may restrict our formahsm to the case ^ = j = 
but aUow, nevertheless, for still arbitrary values of the orbital angular 
momentum £ (which means to consider arbitrary orbital excitations), 
and set (3=1. Then the matrix elements Vij of the power-law potential 
([7]) reduce to 

whereas for the matrix elements Tij of the kinetic energy ([T]) we obtain 

At this point, our primary aim must be to get rid of the hypergeometric 
series F in the above intermediate result. 

• In the ultrarelativistic limit, realized in the case of vanishing mass 
m of the involved particles, that is, for m = 0, the hypergeometric 
series F in may be simplified with the help of the relation 

r{w)r{w -u-v) 

F{U, V] W] 1) = — -— r 

I [w — u) I [w — v) 

for 

7^ 0, -1, -2, . . . , ^{w - u - v) > , 

in order to yield for the kinetic-energy matrix element Too, Eq. (|SD, 
the much more innocent expression 

2[T{i + 2)f 

The resulting upper bounds, iifgo? can be optimized be minimizing 
Hqq with respect to the variational parameter /i. For instance, for 
a linear potential l^(|x|) = a |x|, this minimization procedure thus 
yields 

' (2^ + 3)a 



min//oo = 2r(£+ 2) / / , 

M>o ^ ^ ^^r(^ + |)r(f + | 

In the limit of large orbital angular momenta ^, that is, for ^ ^ 00, 
this minimal upper bound turns out to be not in conflict with the 
experimentally well-established linearity of "Regge trajectories:" 

lim min Hnn = 8 a £ , 

which is in striking accordance with all previous findings [|25|, [26 
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• Fixing the variational parameter /i to the particular value /i = m 
allows us to take advantage from the fact that 

v;w;0) = l , 

whence the kinetic-energy matrix element Too, Eq. (@), reduces to 



Too — 



4^+^ [r(f + 2)]' 



m . 



5.2 Radial Excitations 

On the other hand, considering only states of vanishing orbital angular 
momentum i. e., only states with ^ = 0, confines our investigation to 
the analysis of radial excitations. In this case, we may use the relation 

m 



3 o \ sin[(2w- l)z\ 
V 2 / (2 w — 1} smz 

in order to recast the hypergeometric series F in the momentum-space 
representation '?/^fc,oo(|p|) of our trial states into the form 

^2 



F 



/3 + r + 2 /3 + r 3 



2 '2'p2 + ^2 



(/3 + r + 1) IpI 



X sm 



(/3 + r + 1) arctan — 



Simplifying the momentum-space trial function '?/'fc,oo(|p|) in this way, 

2^ 



V^mo(IpI) = 



k\ 



k 



Vr(2/3 + A; + l) 7r|p 
2Y f k + 2/3 



k — r 



X 



2\ -{f3+r+l)/2 



r(/3 + r + r 



sm (/3 + r + 1) arctan — 



the matrix elements Tij of the kinetic energy immediately become 

4/5+1 



T = 

■J- 1 J 



^r(2/? + z + i)r(2/3 + j + i) TT 



X E E 



r=o s=o r\ s! 



^ + 2/3\ j + 2p 
i-r [ j - s 



X T{(5 + r + l)T{(5 + s + l)Irs , 
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where Irs denotes the only remaining integration, 

3 

dy 



-^rs — 



\ 



X 



cos[(r — s) arctan?/] — cos[(2 + r + s + 2) arctan?/] 

(1 + ^2)(2/3+r+s+2)/2 



This integration may, of course, always be performed by some standard 
numerical integration procedure. However, for /i = m, the integral I, 
simplifies to 



rs 



oo 



Irs = / dy 



cos[(r — s) arctan?/] — cos[(2 /3 + r + s + 2) arctan?/] 

(1 + ^2y2/3+r+s+l)/2 



which, for 2 /3 integer and thus, because of the previous normalizability 
constraint 2 /3 > —1, non-negative, i. e., for the values 2 /3 = 0, 1, 2, . . ., 
may be evaluated with the help of the expansion 



cos(A^ arctan?/) 



1 



N 



N 
n 



cos 



y 



for iV = 0, 1,2, 



with the result 



rs 



1 

2 

X r 
1 

2 

X r 



r 



'2/? + r + s + Ir - si + 1~ 



r 



'2/? + r + S+ Ir 



\r-s\ 

E 

- n 



\r — s\ 
n 



n IT 



cos 



r 2/? + r + s + 



E 

77 = 



n 



'n + 1^ 



r (^2 + r + s + 1 - cos 



n TT 



The case = 0, however, requires special care for the following reason. 
For /? = 0, the integral /qo and therefore also the kinetic-energy matrix 
element Too become singular, as may be read off from the above explicit 
expression for the integral Irs- This singularity may be cancelled by the 
contribution of a Coulomb-type term tz |x|"^ in the power-law potential 
([7]) if the involved coupling constant k takes some particular, "critical" 
value. This cancellation can then be made manifest by observing that 
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oo 



lim 



dy 



1 



cos[(2 + 2/3)arctan7/] ^ J 



V- ' y / '0 

Explicitly, for /3 = 1, the kinetic-energy matrix T = {Tij) is given by 



rj. 128 

1 = m 

IStt 



1 ^ 

7 

V3 11 



\ 



V 



In any case, our approach yields analytic expressions for the matrix 
elements Hij of our semirelativistic Hamiltonian H with an interaction 
potential out of the rather large class given by the power-law form (jT]). 
In principle, the d (real) roots of the characteristic equation (g) may be 
determined algebraically up to and including the case d = 4, entailing, 
of course, analytic expressions of rather rapidly increasing complexity. 
For larger values of the dimension d of our trial space D^, the resulting 
energy matrix, {Hij), may be easily diagonalized numerically, however, 
without the necessity to apply time-consuming integration procedures. 

In order to be able to estimate and appreciate the quality of all the 
upper bounds obtained in this way, we apply the above results to four 
prototype potentials, namely, to 



• the harmonic-oscillator potential, 

y(|x|) =a;|xp 



CO > 



• the Coulomb potential. 



v(W) = 



X 



k: > , 



• the linear potential, 
and 

• the funnel potential. 



y(|x|) = a |x| , a > 



V(|x|) 



+ a|x|, Av>0, a>0. 
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for typical values [j7] of the involved coupling parameters cj, /t, and a. 
The upper bounds on the energy eigenvalues of the lowest-lying radial 
excitations (IS, 2S, 3S, and 4S in usual spectroscopic notation) for the 
harmonic-oscillator, Coulomb, linear, and funnel potentials are shown 
in Tables |I] through g, respectively; the upper bounds on the respective 
energy eigenvalues of just the first orbital excitation (IP again in usual 
spectroscopic notation) for the above potentials are listed in Table |5|. 

Table 1: Energy eigenvalues of the spinless Salpeter equation with harmonic-oscillator 
potential V(|x|) = lu |xp, for the parameter values /i = m = 1 GeV, lu = 0.5 GeV^, 
(3 = 1, and a size dxd of the energy matrix {Hij). Numbers in italics (for small matrix 
sizes) indicate analytically obtained results. All eigenvalues are given in units of GeV. 



State 


1 X 1 


2x2 


25 X 25 


Schrodinger 


IS 


4.2162 


3.9276 


3.8249 


3.8249 


2S 




8.1085 


5.7911 


5.7911 


3S 






7.4829 


7.4823 


48 






9.0215 


9.0075 



Table 2: Energy eigenvalues of the spinless Salpeter equation with Coulomb potential 
y(|x|) = — k/|x|, for the parameter values = m = 1 GeV, /9 = 1, k = 0.456, and 
the size d x d of the energy matrix {Hij). Numbers in italics (for small matrix sizes) 
indicate analytically obtained results. All eigenvalues are given in units of GeV. 



State 1 x 1 2 x 2 25 x 25 



IS 2.2602 2.0539 1.9450 

2S — 3.0702 1.9868 

3S — — 2.0015 

4S — — 2.0238 



For the case of the harmonic-oscillator potential, the corresponding 
Hamiltonian H in its momentum-space representation is equivalent to 
a nonrelativistic Hamiltonian with some effective interaction potential, 
which clearly is reminiscent of that troublesome square-root operator. 
In this form, it is then rather easily accessible to numerical procedures 
for solving a nonrelativistic Schrodinger equation []27[] . For comparison, 
we quote, in Tables |I] and |5|, the eigenvalues obtained along these lines. 

We find a very encouraging, rapid convergence of the upper bounds. 
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Table 3: Energy eigenvalues of the spinless Salpeter equation with the linear potential 
V^(|x|) = a |x|, for the parameter values [0] /i = m = 1 GeV, P = 1, a = 0.211 GeV^, 
and the size d x d oi the energy matrix (Hij). Numbers in italics (for small matrix 
sizes) indicate analytically obtained results. All eigenvalues are given in units of GeV. 



State 1 x 1 2 x 2 20 x 20 



IS 3.0327 2.8034 2.7992 

2S — 4.0767 3.3629 

3S — — 3.8079 

4S — — 4.1905 



Table 4: Energy eigenvalues of the spinless Salpeter equation with the funnel potential 
V^(|x|) = — K/|x|+a |x|, for the parameter values ^ = m = 1 GeV, /3 = 1, k = 0.456, 
a = 0.211 GeV^, and the size dx d of the energy matrix (Hij). Numbers in italics (for 
small matrix sizes) indicate analytically obtained results. All eigenvalues are given in 
units of GeV. 



State 1 x 1 2 x 2 20 x 20 



IS 2.5767 2.5182 2.5162 

2S — 3.4499 3.1570 

3S — — 3.6337 

4S — — 4.0348 



Table 5: Energy eigenvalues for the IP states of the spinless Salpeter equation with the 
harmonic-oscillator potential V^(|x|) = uj |xp, the Coulomb potential V^(|x|) = — k/|x|, 
the linear potential V^(|x|) = a |x|, and the funnel potential V(|x|) = — k/|x| + a |x|, 
respectively, for the parameter values /i = m = 1 GeV, /3 = 1, c<j = 0.5 GeV^, 
K = 0.456, a = 0.211 GeV^, and the size dx d of the energy matrix (Hij). Numbers in 
italics (for small matrix sizes) indicate analytically obtained results. All eigenvalues 
are given in units of GeV. 



Potential 1 x 1 20 x 20 Schrodinger 



Harmonic oscillator 6.5094 4.9015 4.9015 

Coulomb 2.5314 1.9875 — 

Linear 3.2869 3.1414 — 

Funnel 3.0589 2.9816 — 
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6 Summary 

By application of the well-known min-max principle, which represents 
the theoretical foundation of any computation of upper bounds on the 
eigenvalues of self-adjoint operators, to trial spaces spanned by sets of 
basis states which enable us to handle the square-root operator of the 
relativistic kinetic energy T in a satisfactory manner, we demonstrated 
how to derive (for lowest-lying states even analytically!) upper bounds 
on the energy levels of the spinless Salpeter equation with some (linear 
combination of) power-law potentials. Interestingly, in the case of the 
funnel potential, which is the prototype of almost all of the "realistic," 
that is, phenomenologically acceptable, inter-quark potentials used for 
the description of hadrons as bound states of (constituent) quarks, the 
obtained lowest-order approximation to the upper bound on, e. g., the 
ground-state energy is merely some 2 % above the corresponding value. 
Of course, all the bounds derived here may be improved numerically by 
a minimization with respect to the variational parameters introduced. 
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